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Abstract 

In this paper we construct explicitly an infinite number of Hopfions (static, 
soliton solutions with non-zero Hopf topological charges) within the recently 
proposed 3 -|- 1-dimensional, integrable and relativistically invariant field theory. 

Two integers label the family of Hopfions we have found. Their product is 
equal to the Hopf charge which provides a lower bound to the soliton's finite en- 
ergy. The Hopfions are constructed explicitly in terms of the toroidal coordinates 
and shown to have a form of linked closed vortices. 



Recent numerical studies of the Faddeev-Skyrme modified 0(3) sigma- model |]T| 
support existence of static toroidal solutions stabilized by their non-zero Hopf numbers 
0. The emerging string-like structures are quite intriquing and may find applications 
in various physical models of condensed matter physics and gauge field theory. 

In we have introduced the three-dimensional field model which falls into a class 
of higher dimensional integrable models from the point of view of the generalized zero- 
curvature approach 0] . The question posed in was whether this form of integrability 
is linked to the existence of soliton solutions as it is expected from the study of two- 
dimensional integrable models. Our analysis of the model in |^ has indeed revealed 
one non-trivial soliton solution described by a standard Hopf map of unit Hopf index. 
To fully establish connection between integrability and soliton solutions would require 
finding other topological solitons with arbitrary topological charges. This is accom- 
plished in this letter. The equations of motion of the model are solved in toroidal 
coordinates and the space of solutions is found to be represented by a family of maps 
labeled by two integers. The integers count number of times the map winds 
around two independent angular directions. 

The model under consideration is described by the Lagrangean density 

^ = -Vo{Hl^f (1) 

where rjo = ±1, determines the choice of the signature of the Minkowski metric, 
dfiu = Vo diag 15— 1) and the field-tensor H^^ is defined in terms of the three 

component, unit vector field n E as 

H^u = n ■ {d^n x d^n) . (2) 

The action in (|l]) is 0(3) and Poincare invariant. The value 3/4 of the power of if^j, in 
(0) is such that the theory circumvents the usual obstacle of Derrick's scaling argument 
against existence of stable solitons. We are interested in the boundary condition n = 
(0, 0, 1) at spatial infinity. This condition compactifies effectively the Euclidean space 
to the three-sphere S^. Accordingly, n becomes a map: S'^. Due to 7r3(S'^) = 

Z, the field configurations fall into disjoint classes characterized by the value of the 
Hopf invariant Qh- 

Using stereographic projection of 5*^ : 

n = - — I — - (u + u*, -i (u-u*) ,\u 1^ -l) (3) 

1+ I M 1^ ^ ^ 

one obtains 

= T, — < — {d^uduU* - duud^u") (4) 
(1+ I u 1^) 
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where m is a complex scalar field. 
In terms of the vector quantities 



K^ = -{l+\u H^.d^'u = {d^u*d,u) d^u - {d,uf d^u* (5) 



and 



the corresponding equations of motion can be rewritten compactly as 

d'-JC^ = (7) 

The quantity (|^) automatically satisfies the relations 

K^'df.u = Im (K^'df.u*) = (8) 

which play a crucial role in establishing integrability of the model. Indeed, using rela- 
tions (P) and the equations of motion one obtains an infinite number of conserved 
currents given by 

with G being any functional of u and u* only (no derivatives). In reference we 
have analyzed the integrability properties of this theory using the generalized version 
of the zero curvature condition The equations of motion can be represented 
as d^B^ + [A'^ , Bfj_] = 0, with A'^ being a fiat SU(2) connection, and -B^ being an 
operator living in any integer spin representation of SU (2) (for more details see §]). 
The integrability properties emerge due to the fact that 5^ can be put in any integer 
spin representation, and that is a direct consequence of (§) 0, The question we 
address here is whether the above notion of integrability implies existence of the infinite 
many soliton solutions. 

We shall look for time- independent solutions to equations of motion (|^ of the typeQ: 

«(^,e,y') = /(r/)e*("«+"^) (10) 
where we used toroidal coordinates on M'^: 

X = ag~^ sinhr/cosy? , ?/ = ag~^ sinhr/sinyj 

z = ag^^sin^ ; a>0 ; g = cosh — cos ^ (11) 
^We shall take m and n integers in order for u to be single valued. 
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The angles ^, both vary from to 27r and r] varies from to oo. The surfaces of 
constant 7] are toroids that circle the z-axis, ^=constant are spheres and y9=constant 
are half-planes. The corresponding gradient of u becomes in toroidal coordinates: 

Vu = (g/a)e^("^«+"^) (^f{r])e, + imf{r,)e^ + zn/(r/)-^ ^ (12) 

where we introduced the unit vectors spanning the orthogonal toroidal coordinate sys- 
tem having properties ■ Cj = z, j = r/, ip. 

In terms of the vector-field /C = (^K ■ Vm*^ K /{\u\'^ + l) the equations of motion 
(0) become, in the static case, V ■ /C = 0, i.e. /C is solenoidal. Using ( p!OD and (0) we 
obtain for the components of /C: 



sinhi] 



(13) 



/ „2 \-l/4 nl/2fll/2 



with 

(sinhr/)2y + 

where the components of the vector field are defined according to ^ = e^K, + e^V^ + 

Plugging the components of JC into the expression for divergence of the vector field 
in the toroidal coordinates 

^■V^l('^ + '^ + -l-'^-2V,?^ + V,,e^-^-^)) (15) 

a \ or] ot, smhr/ 0(p q smh?; q I 



and using 

drj ^ q ^ dC, ^ q 



dq^ ^ 2 sinh ri dq^ , 2 sin ^ 

Q ; ^ = ^ (16) 



we arrive at: 

d J //' 2m^(sinh?7)^ — cosh?] 

drj (/^ + 1)^ rn?{s\.Vi\iriY + 11? sinh 77 

We take m? > in eq.(|T^). The integration yields: 

1 _ 2/ci coshi] 



+ cosh 7] 

with being the integration constant of the last integration. Imposing boundary 
conditions: 

n — s> (0,0, 1) or I M |— » 00 or / — ^ 00 as (19) 

n — >• (0,0,-1) or M— >-0 or f as r] —>■ 00 (20) 
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one gets 



/ 



2 cosh t] — y Ti? /m? + sinh^ rj 



(21) 



'l + m? /ti? sinh^ rj — cosh 77 

One observes that / depends only on the ratio vr? jr?^ and that > for any value 
of rn? /in?. Note, that taking the limit m — > n in (pT]) and using L'Hospital rule yields: 

1 



lim f = ^ 

sinh ?7 



(22) 



which gives for m = n a solution: 



J.m(i+ip) 

sinh ?7 



(23) 



In the special case of m = = 1 eq.(p3D reproduces the standard soliton solution 
1^, 0. Expression (^) can be written as a composite of the Hopf map together with 
the stereographic map : ^ degree 1. 



A formula: 



E 



:i+ 1 u 



(24) 



describes energy of the static configuration Inserting our solution into ( [2^ ) we 
obtain an expression for the energy E = Em^n of the soliton configuration 



= (27r)^8 • 2^1' r ^ 
Jo 1 



dri sinh 7] ( ^ 

m + 



n 



3/4 



+ /^)^ \ (sinh r7)^y 
which after the ?7-integration yields: 

Em,n = (271)^4 • 2^/^ ^\n\\m\{\n\ + \m\) 



3/2 



(25) 



(26) 



We now turn to calculation of the Hopf numbers. Define functions 
as follows: 



1 
2 



3 
4 



cosh?7 — [r? jvr? + sinh^r] 
(Im/nl — 1) [v? jm? + sinh^ r]' ^ 









COS 1 




sin 1 



1 + (m^/ra^) sinh^ ^ —cosh 77 
(|m/n| — 1) {n^ /m? + sinh^ 77^ ^ 



1/2 



COSn(y9 

— sin 



1,...,4 



(27) 



(28) 
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They provide parametrization of and satisfy : 

^ = ^ = ^ , ^1 +^2=1 (29) 

or equivalently, for n defined from u via relation (§) : 

= Zta.Z ; = (Z*, Z*) ; Z = (^^^ (30) 

Let, furthermore 

Ai='-[z^diZ-d,Z^ Z) (31) 

be a vector potential for the two-form Hij = diAj — djAi. In terms of the toroidal 
coordinates the components of A are 

A^ = qd, ($3$4) ; A^ = -mq ($? + ^t) ; = {^l + '^'4) (32) 

Next, we calculate the vector function B = V ^ A (such that Bi = eijkHjk/2), using 
expression for the curl operator in toroidal coordinates and plugging it into the Hopf 
index defined as : 

QH = ^Jd'xA-B (33) 

for which we now find 



°° 2 
$3 + $f 



-nm (34) 

' ' / 

It is always possible to choose the Hopf index such that it is positive Qh > 0. This 
amounts to the right choice of orientation, which determines the sign of Qh- 

Note, that the following inequahty: 



m| + |?7,| > 2A/|n||m| (35) 



always holds. Hence 



v/|n||m|(|n| + |m|)>y2(|n||m|)=^/^ ^ > {27r'f 4 ■ 2^'/^ \Qh\^/^ (36) 

which agrees with the lower bound result by |^ for the Faddeev-Skyrme model. 

Notice that the solution obtained above implies that the vector field n (see (^)) is 
given by 



ni 



p^l cos (m^ + rnp) ; n2 = -pZ^ sin (m^ + nip) ; n^ = (37) 
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Since / is a function of rj only, one observes that the surfaces of constant are tori. In 
addition, on those surfaces, the hues of constant rii and n2 wind around the tori with 
frequencies in the ^ and (/^-directions given by m and n, respectively. One can check 
that ^3 falls monotonically from 77,3 = 1 at r/ = to 723 = — 1 at r/ = cxd. In addition, 
the bigger the ratio m/n is, the faster if performs that flip. Therefore the size of our 
soliton decreases with the increase of m/n. 

The Hopf index of our solution can alternatively be calculated in the following way. 
The solution provides a mapping of the spatial into defined by Yli=i '^'^ = 1- 
Then, (|^) and provide the Hopf map S'^ S^. The Hopf index is given by 
the linking number of the pre-images, under the Hopf map, of any two points of S^. 
Consider the point n = (0, 0, —1) which corresponds in to 1] ^ 00. For this value 
of rj, $i,$2 go to zero and $3 = cosn(f and $4 = —sinn(f. On we find the 
circle $1 + $4 = cos^ mp + sin^ rup = 1 of radius 1 and as ip varies between and 
27r the preimage wraps | n \ times around the ip direction. Similarly the preimage of 
n = (0, 0, 1) corresponds to rj = 0. For this value of rj, $3, $4 go to zero and $1 = cosm^ 
and $2 = sinm^. On 5*^ we find the circle $^ + $2 = cos^ + sin^ = 1 of radius 
1 and as ^ varies between and 271 the preimage wraps | m | times around the ^ 
direction. Since these two circles intersect the linking number is | nm \ . That is indeed 
the Hopf number calculated above. 

Finally, let us mention that other models circumventing Derrick's theorem have 
been proposed previously in the literature @, ^. A common feature is that, like the 
model considered in the present paper, the Lagragians are non-polynomial functions 
of the fields and their derivatives. The corresponding solutions have been constructed 
and in the special case of a soliton with the unit Hopf charge was obtained. To our 
knowledge, an infinite number of soliton solutions to the field theoretical equations of 
motion with the Hopf charges bigger than unity were not obtained previously in the 
literature in an exact form. 
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